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Abstract 



Various astronomical observations have been consistently making a strong 
case for the existence of a component of dark energy with negative pressure 
in the universe. It is now necessary to take the dark energy component (s) 
into account in gravitational lensing statistics and other cosmological tests. 
By using the comoving distance we derive analytic but simple expressions for 



^ ' the optical depth of multiple image, the expected value of image separation 



and the probability distribution of image separation caused by an assemble 
of singular isothermal spheres in general FRW cosmological models with dark 
energy component(s). We also present the kinematical and dynamical proper- 
ties of these kinds of cosmological models and calculate the age of the universe 
and the distance measures, which are often used in classical cosmological tests. 
In some cases we are able to give formulae that are simpler than those found 
elsewhere in the literature, which could make the cosmological tests for dark 
energy component (s) more convenient. 
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I. INTRODUCTION 

The Hubble expansion, the Cosmic Microwave Background Radiation(CMBR) and the 
primordial Big Bang Nucleosynthesis(BBN) are three cornerstones of the standard hot Big 
Bang cosmological model. Recently these three kinds of observations have been making a 
strong case for the existence of a nearly uniform component of dark energy with negative 
pressure: More than twenty experiments of CMBR measurements have now covered three 
orders of magnitude in multipole number, and are beginning to define the position of the first 
acoustic peak at a value that is consistent with a flat universe [|1|; However, the deuterium 
abundance measured in four high redshift hydrogen clouds seen in absorption against distant 
quasars combined with baryon fraction in galaxy clusters from X-ray data ||^ give a low 
matter density universe, Qm ~ 0.3 — 0.4; The discrepancy between the low matter density 
and a flat universe is observationally resolved by high redshift la type supernovae (SNela) 
observations @], which implies an accelerating universe driven by a dark, exotic form of 
energy component. For a recent review of observational evidence for the dark energy com- 
ponent, reader is referred to Ref. |^. It seems that determining the amount and nature of 
the dark energy is emerging as one of the most important challenges in cosmology. 

The simplest possibility for the dark energy component is cosmological constant A or 
equivalently vacuum energy 0. In the past, a nonzero cosmological constant has been 
advocated and then with the improved observational data discarded several times [0]. Due 
to this checked history and the difficulty in understanding the observed A in the framework of 
modern quantum field theory, now most physicists and astronomers prefer other candidates 
for dark energy, including a frustrated network of topological defects such as strings or walls 
P] and an evolving scalar field referred to by some as quintessence P] etc. As shown in 
literatures, it is difficult to discriminate against these different possibilities either by the 



SNela data alone [T^ or only by the CMBR data [TT|. This led some authors to consider 



the combination of the SNela measurements with the anisotropy of CMBR [|T2| or the large 



scale structures 13 



Because of the potential importance of probing the dark energy component and uncer- 
tainties of present constraints, it is worthy of examining various cosmological tests for dark 
energy component (s). It has been long shown that gravitational lensing statistics is an effi- 



cient tool for determining cosmological parameters |]14|,|15| . Some authors have even given the 
general expressions for the optical depth and average image separation in general Friedman- 
Robertson- Walker (FRW) cosmological models ||T6[. But these expressions are complicated 
and thereby hard to apply in practice. Furthermore, they didn't involve general dark energy 
component (s). We reconsider here gravitational lensing statistical properties, such as the 
optical depth of multiple image, the expectation value of image separation and the prob- 
ability distribution of image separation, in general FRW cosmologies with various matter 
and energy components. We also calculate the age of the universe and distance measures, 
such as the luminosity distance and the angular diameter distance, which are often involved 
in classical cosmological tests. In some cases we are able to give formulae that are simpler 
than those found elsewhere in the literature, which could make the cosmological tests for 
dark energy component (s) more convenient. 

The paper is organised as follows: In section 2, we summarize kinematics and dynamics of 
general FRW cosmologies with dark energy component (s); In section 3, we derive expressions 
for the age of the universe and the various distance measures; In section 4, using the comoving 
distance, we generalize and simplify lensing formulae; Results and discussion are given in 
section 5. 

II. KINEMATICS AND DYNAMICS OF GENERAL FRW COSMOLOGIES WITH 

DARK ENERGY COMPONENT(S) 

For general FRW cosmologies, the metric of spacetime is described by (in the c = 1 unit): 

ds^ = -de + i?2(t) [dx^ + f{x){de^ + sin' ed(f)^)] , (1) 

where f{x) = X for ^ A^t universe {k = 0), f{x) = sinx for a closed universe {k = +1), 
and f(-v) = sinh v for an onen universe (k = —1). Defining the scale factor ait) = R(t)/Rn. 



a = 1 today, and the Friedmann equation takes the form 

fa\'^ _ SttGx^ „ k 



(2) 



-^Ei{p^ + Spi) 



3 

where the dot represents derivatives with respect to t, and i includes all components of 
matter or energy in the universe, e.g., i = m for the total nonrelativistic matter, r for the 
total radiation component, A for the cosmological constant (pa = A/(87rG')), c for the cold 
dark matter, h for the hot dark matter, b for the baryons and x for some unknown exotic 
component. If the effective equation of state for the i-th component is parameterized as 
uJi = Pi/ pi, its density scales as pi oc a~"' where Ui = 3(1 +uji). For instance, nonrelativistic 
matter scales as pm oc a^^ while relativistic matter, such as radiation, changes as pr oc a~^, 
and vacuum energy(cosmological constant) is invariant {p\ oc a°) as the universe expands. 
Eq.(^ tells us that all components of matter or energy and a curvature term conspire to 
drive the universal expansion. It is convenient to assign symbols to their respective fractional 
contributions at the present epoch. Defining parameters 

8ttG n _ -k 

where Hq is the Hubble constant and zero subscripts refer to the present epoch. Eq.(^ 
implies a relation 

l = ^a + fifc. (4) 

i 

The parameters fij, Vt^ can be used to rewrite equation (^ at a general time 

i|-^o(i: ".«---' + ^,.-)"^ (5) 



III. THE AGE OF THE UNIVERSE AND DISTANCE MEASURES 

Equation (|D can be easily used to calculate the age of the universe. The age is the 
inteErral of the time un to the nresent enoch in terms of the normalized scale factor a: 



Hoto 



da 



(6) 



By a trivial change of the integral variable from a to redshift z, the age expression reduces 



to 



Hnt, 



O'-O 



dz 



[l+z)[j:rn,{l + zY(^+'^^) + Qk{l + 



N2lV2 



(7) 



It has been shown that the natural cosmological distance for the analysis of gravitational 
lensing statistics is the conioving distance |17|. The distance a light ray travels can be 
calculated as follows. Light rays follow null geodesies where ds"^ = so that dt"^ = R^a?dx^. 
With Eq.dID and the relation a{t) = (1 + z)~^, the comoving distance is 



v/E7mi+^ 



1+^,) 



x = \ 






(8) 



la 



1 1/2 



/o 



dz 



{k = ±l). 



Most of classical cosmological tests involve the standard candle (e.g. SNela) or standard 
ruler (e.g. radio source). Therefore both luminosity distance D^ and angular diameter 



distance D are often used. They are related to the comoving distance as |jT8 



D^^Rlf{x)lR{t). D^^R{t)f{x). (9) 

They are simply related to the measurable, the redshift z, through Eq.(|^) with the following 
form 

H^\l + z)x k = 0, 

k = +1, (10) 



D' 



'l+zYD 



2nA 



1-1/2, 



Ho\l + z)\nk\~'^^^smhx k = -l 



Finally, one often needs to calculate volume element of the universe. In order to take 
into account the cosmic expansion, it is more convenient to use comoving volume than the 
traditional physical volume. Within the shell dx at x? dV reads 

X^dx k = 0, 



dV = A'KRlf\x)dx = 47riJo" 






' sin^ xdx k = +1, 
' sinh xdx k = —1. 



(11) 



IV. GRAVITATIONAL LENSING STATISTICAL PROPERTIES OF SINGULAR 

ISOTHERMAL SPHERES 

A. Optical depth for multiple image 



First of all, we consider the lensing cross-section [IS] due to a specific galaxy. Following 



Ref. [IS], we model the mass density profile of the total galaxy matter as the singular 



isothermal sphere (SIS). The dimensionless cross-section of multiple image for a point source 



located at Zs produced by a single SIS galaxy at Zd is [20 



a = nel, OE^ina'^, (12) 

where D^ and D^^ are the angular diameter distances from the observer to the source and 
from the lens to the source respectively, 6e is the angular radius of Einstein ring and a is 
the velocity dispersion of the lensing galaxy. Using the comoving distance gives 



a = levrV^ 



fiXs -XdY 



(13) 



f{Xs) 

Now, let's consider the contributions of an ensemble of galaxies having different lumi- 
nosities and redshifts. The present-day galaxy luminosity function can be described by the 



Schechter function 21 



UL)dL = 0*(L/L*)-"' expi-L/L*)diL/L*), (14) 

where i indicates the morphological type of galaxies: i=(E, SO, S). The above expression can 
be converted into the velocity dispersion distribution through the empirical formula between 
the luminosity and the central dispersion of local galaxies L/L* = (cr/cr*)^\ Finally, the 
optical depth of multiple image by galaxies at redshifts ranging from to Zg for the distant 
sources like quasars at Zs is 

rizs)= ( E fAt{zs), (15) 

\i=E,SO,S / 

The parameter Fi represents the effectiveness of the i-th morphological type of galaxies in 
nroducine double images 111 Oil. which reads 



Fi = 16n^Ho-'{no,a^) = IGTr^i/o^'^hi (ha*)^ J {1/1,^+^/^^ exp{L/L,)dL/L, , (16) 

where 7^ is the galaxy morphological composition and bi is the velocity bias between the 
velocity dispersion of stars and of dark matter particles [^. The above equation can be 
further written as 



F, = lQn'Ho-'(l)Hihcr:rr{-a, + A/g, + l), if L G (0, 00) 



(17) 



if the integral is performed from to 00. In practice, the galaxy luminosities have the 
minimum and maximum limits, and, therefore, Eq.(p!7D is the maximum estimate of F,. The 
Zs dependent factor T{zs) is 



T(z.. 



(HoRq) 



Xs 



fiXs - Xd) 

fiXs) 



f{Xd)dXd- 



(18) 



For general FRW cosmologies with various matter and energy components, an analytic 



expression is found |17 



Tiz.. 



2G 
30' 



10 r^/^ 



10 r3/2 



(A; = 0), 
{k = -l) 



(19) 



L*,oti]cr*,gi]^i)oiga\ax.- 



\{l + ?,coi'^Xs)Xs - IcoiXs 
|(-l + 3coth^Xs)Xs- fcothx. 

where Xs can be calculated through Eq.(||). 

As pointed out by Kochaneck [^, although the parameters 
ies are updated due to recent surveys [^ there are still large uncertainties with them. 
Therefore we treat F = J2i^i as normalized factor. One may expect to determine both the 
dark energy amount and its equation of state. It is necessary to consider general flat universe 
models [flm] ^x, ^x) with Qm + fi^; = 1. Generally, the optical depth of gravitational lensing 
depends on the amount of dark energy Q^ as well as its equation of state Ux- The larger 
the dark energy amount is, the higher the gravitational lensing probability will be; the more 
negative the dark energy pressure is, the higher the optical depth will be. 



B. Average image separation 

In addition to optical depth, other interesting quantities for lensing statistics are the 
mean image separation produced by a specified SIS at a given z^ with velocity dispersion 
a and its average over the distributions of lens velocity dispersion and lens redshift. The 
angular separation of the two images produced by SIS is 

A^ = 2eE = 87ra'^^%^^ (20) 

fiXs) 

which doesn't depend on impact parameter, and hence implies 

M = A^ for SIS (21) 

where A^ is the mean image separation for the SIS lens at Zd with velocity dispersion a. 
The average image separation is obtained by 

<A9>=-fAedT (22) 

We take into account the lens redshift distribution as well as the lens velocity dispersion 
distribution. The latter was generally neglected in literatures. Then the average image 
separation reduces to 

<M>=M Y. ^'|0(^s), (23) 

^ \i=E,SO,S J 

where 

F^ = 1287r^i/o"'K.a6) = l28n*Ho-Y^l^ {half j [L/ UY^^''/^^ exp{L/L,)dL/L, , (24) 
if the above integral is performed from to oo, it can be further written as 

Fl = 1287r^/7o"'0*7i {k(r*f r(-a + 6/gi + 1), if L G (0, oo) , (25) 

and the Zg dependent factor Q{zs) is 



e(^,) = {HoRof r 

.In 



fixs-XdV 



f{Xd)dXd. (26) 



For general FRW cosmologies with various matter and energy components, an analytic 
expression is found: 



e{zs 



2G 
60' 



\n,r'/'^^sm-^xs 



{A-3sin^Xs) +cosXs(-4 + sin2xs 



\n, 



-3/2 2 



Yg sinli Xs (-4-3 sinli Xs) + cosh Xs(4 + sinh Xs 



{k = +l] 
{k = -l] 



(27) 



where Xs can be calculated through Eq.(||) as before. If we define A^* = Y,i F'll Si Fi^ which 
is determined completely by the intrinsic characteristic parameters of galaxies describing its 
statistical properties, the average image separation then becomes 



<M>=M,e{zs)/T{zs 



Ae, 



yi=E,SO,S 



.i=E,SO,S 



(28) 



Thus for an open or closed universe, the average image separation for the SIS lens depends 
on Xs (or the source redshift), while for a flat universe, it doesn't but takes a constant value 
< A9 >= A9^/2. This fact has been suggested to test directly the curvature of the universe 



24 . 



C. Probability distribution of image separation 



Another useful quantity of lensing statistics is the probability distribution of image sep- 
aration, which is reconsidered here for general cosmologies with various matter and energy 
components. The differential optical depth for multiple imaging a background object like a 
quasar by lenses at position between Xd and Xd + dXd and with luminosity from L to L + dL 



IS 



d'n = -fiML)dL a Rlf\xd)dxd. 



(29) 



By change of the variable from luminosity to angular separation, the differential optical 
depth becomes 



(30) 



f{Xs) 



9./2 



■exp -[87r(6,a„)r^'/' ^%^ ^(^''"\ ^Hfvf^ RlPiXa 



fiXs-Xd) 



fiXs) 



gi{ai~l) 



where we have taken into account the angular resolution bias by the function S{A6), which 



varies for different lens observations and surveys ^^. Then, the probability distribution of 
image separations becomes 

which in general case depends on Xs (or Zg, the redshift of the source). However for a 
flat universe, the above integral scales as oc Xs- Combining Eq.(|T^) and |TP|, one concludes 
that the probability distribution of image separations is independent of Xs, regardless of the 
functional form of the angular resolution bias |2^. Hence the cumulative distribution of 
image separations can be easily calculated from dP/dA6 as 

fAe [ (jP \ 



V. RESULTS AND DISCUSSIONS 

As shown above, it is easy to analytically calculate the optical depth for multiple image 
(Eq.([l5|) and (|l^)), average image separation (Eq.f^SD, (|T^) and (|3)) and probability dis- 
tribution of image separations (Eq. (|3l|) ), as well as the age of the universe (Eq.(^) and the 
luminosity (angular diameter) distance (Eq.([lOD), if the matter and energy components in 
the universe are given in terms of [Vti^uoi). It is shown that the lensing statistical properties 
depend sensitively on the dark energy amount and its equation of state and hence provide 
an independent probe for the dark energy. 

In the above lensing calculations, we have assumed that the comoving number density 
of galaxies is constant. However, this may not hold true for the realistic situation. The 
influence of ffalaxv evolution on the lensinff statistics should also be taken into account 1 12 711. 



We have considered this effect by using the galaxy merging model proposed by Broadhurst et 
al. (1992), since the scenario of galaxy merging can account for both the redshift distribution 
and the number counts of galaxies at optical and near-infrared wavelengths ||2^ . There are 
two effects arising from the galaxy merging: The first is that there are more galaxies and 
hence more lenses in the past. The second is that galaxies are typically less massive in the 
past and hence less efficient as lenses. As a result of two effects, the total cross-section 
remains roughly unchanged [^ . 



Although the evidence for the existence of a nearly uniform component of dark energy 
with negative pressure in the universe has been consistently provided by several astronomical 
observations, we at present know a little about it. It needs more cosmologcial constraints 
on dark energy to determine its amount and nature. The simple formulae we derived here 
could make the cosmological tests for dark energy component (s) more convenient. 

ACKNOWLEDGMENTS 

I thank X. P. Wu for careful reading of the manuscript and helpful discussions. This 
work was supported by the National Natural Science Foundation of China, under Grant No. 
19903002. 



REFERENCES 

[1] C. Lineweaver, Astrophys. J., 505 (1998) L69. 

[2] S. Buries, D. Tytler, Astrophys. J., 499 (1998) 699; S. Buries, D. Tytler, Astrophys. J., 
507 (1998) 732. 

[3] U.G. Briel, J.P. Henry, H. Bohringer, Astron. Astrophys., 259 (1992) L31; S.D.M. White 
et al.. Nature, 366 (1993) 429; J.J. Mohr, B. Mathiesen, A.E. Evrard, Astrophys. J., 
517 (1999) 627. 

[4] S. Perlmutter et al., Astrophys. J., 517 (1999) 565; A.G. Riess et al., Astron. J., 116 
(1998) 1009. 



[5] M.S. Turner, |astro-ph/98Tl454 . 



[6] S. Weinberg, Rev. Mod. Phys. 61 (1989) 1; S. Carroll et al., Ann. Rev. Astron. Astro- 
phys., 30 (1992) 499; L.M. Krauss, M.S. Turner, Gen. Rel. Grav., 27 (1995) 1137; J.P. 
Ostriker, P.J. Steinhardt, Nature, 377 (1995) 600. 



[7] J.D. Cohn, Living with Lambda, |astro-ph/980712^ ; V. Sahni, A. Starobinsky, Int. J. 



Mod. Phys. D, in press ( |astro-ph/9904398|) . 



[8] A. Vilenkin, Phys. Rev. Lett. 53 (1984) 1016; R.L. Davis, Phys. Rev. D35 (1987) 3705; 
M. Kamionkowski, N. Toumbas, Phys. Rev. Lett. 77 (1996) 587. 

[9] B. Ratra and P.J.E. Peebles, Phys. Rev. D37 (1988) 3406; J. Frieman et al, Phys. Rev. 
Lett. 75 (1995) 2077; K. Coble et al, Phys. Rev. D55 (1995) 1851; R. Caldwell et al, 
Phys. Rev. Lett. 80 (1998) 1582. 

[10] P. Garnavich et al, Astrophys. J., 509 (1998) 74. 

[11] G. Huey et al., Phys. Rev. D59 (1999) 063005. 

[12] M.S. Turner, M. White, Phys. Rev. D56 (1997) R4439. 



[13] S. Perlmutter, M.S. Turner, M. White, Phys. Rev. Lett., 83 (1999) 670. 

[14] E.L. Turner, Astrophys. J., 365 (1990) L43; M. Fukugita et al., MNRAS, 246 (1990) 24; 
L.M. Krauss, M. White, Astrophys. J., 394 (1992) 385; C.S. Kochanek, Astrophys. J., 
466 (1996) 638; A.R. Cooray, Astron. Astrophys., 341 (1999) 653; R. Quast, P.Helbig, 
Astron. Astrophys., 344 (1999) 721; M. Chiba, Y. Yoshii, Astrophys. J., 510 (1999) 42; 
I. Waga, A.P.M.R. Miceh, Phys.Rev.D59 (1999) 103507; Y. Melher, Ann. Rev. Astron. 
Astrophys., 37 (1999) 127; Z.H. Zhu, Mod. Phys. Lett. A15 (2000) 1023; Y.-C.N. Cheng, 
L.M. Krauss, Int.J.Mod.Phys.Al5 (2000) 6971. 

[15] X.P. Wu, F. Hammer, MNRAS, 262 (1993) 187; X.P. Wu, S. Mao, Astrophys. J., 463 
(1996) 404; T. Hamana, T. Futamase, MNRAS, 286 (1997) L7; M. Bartelmann et al., 
Astron. Astrophys., 330 (1998) 1. 

[16] J.R. Gott, M.G. Park, H.M. Lee, Astrophys. J., 338 (1989) 1; C.S. Kochanek, MNRAS, 
261 (1993) 453. 

[17] Z.H. Zhu, Astron. Astrophys., 338 (1998) 777; Z.H. Zhu, Mod. Phys. Lett. A15 (2000) 
1023. 

[18] S. Carroll et al., Ann. Rev. Astron. Astrophys., 30 (1992) 499. 

[19] E.L. Turner, J.P. Ostriker, J.R. Gott, Astrophys. J., 284 (1984) 1. 

[20] Schneider P., Ehlers J., Falco E. E., 1992, Gravitational Lenses (Springer Verlag, Berlin); 
X.P. Wu, Fund. Cos. Phys. 17 (1996) 1. 

[21] Peebles P. J. E., 1993, Physical Cosmology (Princeton Univ. Press, Princeton). 

[22] C.S. Kochanek, Astrophys. J., 466 (1996) 638. 

[23] J. Loveday, B.A. Peterson, G. Efstathiou, S.J. Maddox, Astrophys. J., 390 (1992) 338; 
R.O. Marzke, M.J. Geller, J.P. Huchra, H.G. Corvin, Astro. J., 108 (1994) 437. 

[24] M.G. Park, J.R. Gott, ApJ 489 (1997) 476; P. Helbig, MNRAS 298 (1998) 395. 



[25] A.O. Jaunsen, M. Jablonski, B.R. Pettersen, R. Stabell, Astron. Astrophys., 300 (1995) 
323. 



[26] J.R. Gott, M.G. Park, H.M. Lee, Astrophys. J., 338 (1989) 1; M.G. Park, 1996, |iitm 



ph/9609204 . 



[27] S. Mao, C.S. Kochanek, MNRAS, 268 (1994) 569; D. Jain, N. Panchapakesan, S. Ma- 
hajan, V.B. Bhatia, Mod.Phys.Lett.Al5 (2000) 41. 

[28] T. Broadhurst et ah. Nature, 355 (1992) 55. 

[29] Z.H. Zhu, X.P. Wu, Astron. Astrophys., 326 (1997) L9. 



